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Abstract Magnetic Field extrapolation is an alterna- 
tive method to study chromospheric and coronal mag- 
netic fields. In this p aper, two sern i -analy tical solu- 
tions of force- free fields (|Low and Loul . 119901 ) have been 
used to study the errors of nonlinear force-free (NLFF) 
fields based on force-free factor a. Three NLFF fields 
are ex trapolated b y app roximate vertical integration 
fAVl HSong et all (l2006l). bou ndarv integral equation 
(BIE) Yan and Sakurai ( 200Gf ) and optimization (Opt.) 
IWiegelmann (2004) methods. Compared with the first 
semi-analytical field, it is found that the mean values of 
absolute relative standard deviations (RSD) of a along 
field lines are about 0.96-1.19, 0.63-1.07 and 0.43-0.72 
for AVI, BIE and Opt. fields, respectively. While for 
the second semi-analytical field, they are about 0.80- 
1.02, 0.67-1.34 and 0.33-0.55 for AVI, BIE and Opt. 
fields, respectively. As for the analytical field, the cal- 
culation error of {\RSD\) is about 0.1 ^ 0.2. It is also 
found that RSD docs not apparently depend on the 
length of field line. These provide the basic estima- 
tion on the deviation of extrapolated field obtained by 
proposed methods from the real force-free field. 
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1 Introduction 

The magnetic field plays key roles in a variety of dynam- 
ical processes, particularly in eruptive phenomena such 
as filament eruptions, flares and coronal mass ejections. 
The topological structure is an important properties of 
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spatial magnetic field. Therefore it is important to un- 
derstand three-dimensional properties of magnetic field 
in order to elaborate efficient theoretical models, which 
will contribute to understand solar activity. At present, 
due to the restrictions of observational technique, we 
can not get the accurate information on magnetic field 
in the hot and tenuous corona. This can only be done 
at the level of the cooler and denser photosphere. In 
spite of that, the extrapolation of magnetograms, mea- 
sured maps of the photospheric or low chromospheric 
field, into the corona is the most important technique to 
supply the missing information about the spatial mag- 
netic field. The field extrap olation is based on force- 
free assumption ( Alvl . Il989l ). which assumes that the 
corona is free of Lorentz forces. For the low-/? corona 
where the plasma is tenuous, the force-free assumption 
is appropriate. The force-free field obeys the following 
equations: 



V X B = a{r)B, 



VB = 0, 



(1) 



(2) 



where a is called force-free factor. If a = 0, the equa- 
tions represent a potential field (a current- free field). If 
a = constant, they describe a current-carrying linear 
force-free (LFF) field, and if a = /(r) it is a general 
NLFF field. This two equations also imply that: 



B • Va = 0. 



(3) 



Equation ([3]) demonstrates that a is invariant along 
field lines of B . The scalar a is in general a function of 
space and identifies how much current fiows along each 
field line. Thus, the deviations of a along one field line 
at some extent can demonstrate performance of field 
extrapolation. 
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At present, the potential (a — 0) and linear force- 
free {a = constant) field extrapolations reached a 
mature development, but they describe the magnetic 
field above the photosphere in a very restricted man- 
ner. According to the practice conditions of so- 
lar magnetic field, it is not a current-free field and 
NLFF should be more reasonable. Recently, sev- 



B„ 



^2ix,y,z)F2ix,y,z), 



= ^3{x,y,z)F3{x,y,z), 



(5) 



(6) 



posed (e.Q., Wu et all 


1990; Cuoerman et al. 




1990; 


Demoulin et ai. 1992; 


Mikic and McClvmont 


1994; 


Roumeliotis. 1996: Amari et ai. 1997; Sakura] 




1981; 


Chodura and Schlueter, 1981; Yan and Sakura 




2000; 


Wheatland et ai. 200(i WieKelmann,i2004;,SonE et ai. 


20061 He and Wansl 20081). The performance of above 



extrapolatioi i methods have been studied in s everal 
papers (e.g.. ISchriiyer et al. . 20061 Ainari et ai . 2006 



assuming the second-order continuous partial deriva- 
tives in a certain height range, 0< z <H (H is the 
calculated height from the photospheric surface). In 
Equations (H])-®, £,1,^2 and ^3 mainly depend on z and 
slowly vary with x and y, while Fi,F2 and F^ mainly 
depend on x and y and weakly vary with z. After con- 
structing the magnetic field, the following integration 
equations. 



2006 Valori et al 
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2011 



Regnier a/j.l2004tlYan and Lil.l2006l : |Wiegelmann adl — i = — 



and 



Liu et al\ 1201 ll ) , compar 



■ l2007tlRegnier et a4f2007tlDeRosa et a J |'^^ 



ing to models or observations (e.g., the appearance of 
extreme ultraviolet (EUV) and X-ray or loops) those 
methods above mentioned can give reasonable results 
at the level of macroscopic structure. 

Since the field extrapolation is an important tech- 
nique to study the coronal magnetic field, the per- 
formance of field extrapolation become special impor- 
tance. From Equation ([3]), it can been seen that though 
a may change from one field line to another, it must be 
a constant along one field line. Knowing a is an impor- 
tant parameter worthy to be researched, in this paper 
the distributions of a between extrapolated field and 
semi-analytical field are compared; the calculation er- 
ror of a is discussed; the performance of a along field 
line are studied specially. 

The paper is organized as follows: firstly, the de- 
scription of extrapolation methods and semi-analytical 
field will be introduced in Section [2] secondly, the re- 
sults of distribution comparison and calculation error 
of force-free factor a are shown in Section |3| at last, 
the short discussions and conclusions will be given in 
section 



2 Extrapolation methods and Semi-analytical 
field 

2.1 Approximate vertical integration method 

The approximate vertical integration (AVI) method 
I Song et ali . 120061 ) w as improved from the vertical inte- 



gration proposed by Wu et al. ( 1990l ). In this method 



it is assumed that the magnetic field components is 
given by the following formula. 



dz 



dy 

dB. 



+ aBy, 

- aB^, 
dB, 



(7) 

(8) 
(9) 

dx dy ' ^^^^ 

can be used to calculate extrapolated field. There ex- 
ists singularity problem due to differential operation for 
AVI method, h ence we sn i ooth d ata locally where exists 
singularity ( cf. iLiu et ali . 120111 ) . 



d2 



dx 



dy ' 



aB, = 



dB„ dB,, 



2.2 Boundary integral equation method 

The bou ndary integral eq u ation (BIE) method pro- 
posed by lYan and Sakurail (l2000l) . which uses the in- 



tegration function to extrapolate the magnetic field. In 
this method, an optimized parameter A, which is the 
function of spatial position x, must be found through 
iteration. The integral 



B(a;i,2/i,Zi) 



Zi[\rs\n[\r) -|- cos(Ar)]Bo(x, y, 0) 
27r[(x - x,Y + [y- y^f + z^Y'^ ' 



(11) 



is used to calculate the magnetic field, where r = [{x — 
+ (y ~ Vi)^ + -2^1]^^^ and Bo is the magnetic field of 
photospheric surface. This method is to find suitable 
values of A through iteration and try to make sure that 
the extrapolated f i eld is for ce-free and divergence-fre e 
(cf. iHe and Wanei bood and lli. Yan. and Songl . 12004) . 



B, 



£i{x,y,z)Fi{x,y,z), 



(4) 
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2.3 Optimization method 



The optimization (Opt.) method proposed b y l Wheatland 
iiQOO) and developed by I Wiegelmam j (j2004l ) consists in 
minimizing a joint measure for the normalized Lorentz 
force and the divergence of the field, given by the func- 
tion, 



field lines of extrapolated fields (AVI, BIE or Opt. field) 
are basically consistent with that of SAFl, and the field 
eaiafcj nf AVI and BIE fields are very similar. Hence the 



L = [ a;(x,y,z)[B-2|(VxBxB)|2 + |V 
Jv 



■B\^]d^x, (12) 



where Ltj{x,y,z) is a weighting function related posi- 
tion. It is clear that (for w > 0) the force-free equa- 
tions are fulfilled when L is equal to zero. This method 
involves minimizing L by optimizing the solution func- 
tion B(x,t) through states that are increasingly force- 
and divergence-free, where t is an artificial time-like pa- 
rameter introduced. 

2.4 Semi-analytical NLFF Field 



Low and Lou (Il990l) have presented a class of axis- 
symmetric NLFF fields in the spherical coordinate sys- 
tem. These NLFF field can be used to test the perfor- 
mance of field extrapolation by shifted under two steps 
of Cartesian coordinate system transformation, rotat- 
ing by an angle (j) around the y-axis and moving the 
origin to a distance I along the z-axis. A special class 
of NLFF fields denoted B(n, m) can be obtained in the 
spherical coordinate system. For different n and m, it 
can give different distributions of magnetic field, which 
meet the requirements of divergence-free and force-free 
equations. Then we should specify two parameters / 
and 4> for coordinate transformation. In this paper we 
choose such two classic NLFF fields as semi-analytical 
fields: 

SAFl: the semi-analytical field with n = 1,™ ~ 1, 
I = 0.3, and ^, set a; e [-0.5,0.5], y G [-0.5,0.5] 
and z G [0, 1] in the Cartesian coordinate system. 

SAF2: the semi-analytical field with n = 3,to = 1, 
47r 

I = 0.3, and (/)=—, set X G [-0.5,0.5], y G [-0.5,0.5] 
5 

and z G [0,1] in the Cartesian coordinate system. 

The mesh is 64 x 64 x 64 for those two fields. Be- 
cause a will be studied mainly, the values of analyt- 
ical a (aAna. ) of S AFl and SAF2 are calculated (c/. 
Low and Loul . 1990 ) and saved in 3D position. 



3 Results 
3.1 SAFl 

Figure ID shows the field lines of SAFl, AVI ,BIE and 
Opt. fields. It can be found that the distributions of 



structure of spatial field obtained from field extrapo- 
lation can present for that of semi-analytical field at 
some extent. In order to evaluate the performance of 
extrapolated field, we calculate the average ((|-F|)) and 
maximum (Fmax) Lorentz- force (J x B//io, fJ-o is per- 
meability) and aj, d efined by equation ([T3|) same as 
Schriiver et all ( 20061 ). to test the extent of force- free 



for each extrapolated field. The results of these phys- 
ical quantity are given in table 1 for SAFl and the 
corresponding extrapolated fields, here the pixel size is 
assumed 1 arcsec. It can be found that angles between 
J and B of Opt. field (comparable to SAFl at some 
extent) are smaller than those of AVI and BIE fields. 



E 



|J, X B, 



(13) 



Because a 3D field is available, force-free factor a 
can be calculated individually from the formulas ((T4| - 
(fT6|) . which are obtained by expanding Equation 
Because a is the scalar depended on spatial position, 
ai, a2 and as should be the same for a given posi- 
tion. Figure [5] shows the images of ai, a2 and a^, mag- 
netic field components and current components ( J^., Jy, 
Jz) labeled for SAFl at z = 0, where the red lines on 
each grey-scale map of magnetic field component are 
the neutral line of each component labeled. Because a 
is calculated from the semi-analytical field, this figure 
demonstrates that there will be unavoidable calculation 
errors when a is calculated from magnetic field. The 
evident calculation errors are most likely located near 
where the magnetic field components reversed, which 
are consistent among ai, and ^3. From the formu- 
las P^ - p6p . it can be seen that there are singularity 
problem when magnetic field component approaches to 
zero, hence the evident calculation errors correspond to 
the neutral line of magnetic field components (red lines 
labeled) . Since a is a scalar and the function of spatial 



Table 1 The results of the average ((j-Fj)) and maximum 
{Fmax) Lorentz-force and dj for SAFl and the correspond- 
ing extrapolated fields. 





m 

{G'^/M)x 10-" 


{G^/M)x 10-12 




SAFl 


2.2 


0.8 


0.14 


AVI 


5.6 


9.6 


0.72 


BIE 


3.7 


5.4 


0.45 


Opt. 


2.3 


2.4 


0.21 
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position, it is reasonable to combine ai, and 0:3 to 
minimize the calculation errors. 



Oil 



'^y/By, Jy 





dBy 


dy 


dz 


dBx 


OB, 


dz 


dx 


dBy 


dBx 



Bx, 



B,, 



(14) 



(15) 



(16) 



In order to compare a, the combined value from ai, 
a2 and 03 is also studied in this work. The method 
to realize this is that: to get the average of two val- 
ues of a {ai, a2 or a3) whose values are closed, which 
means that for a given position a deduced from one 
formula ([Ml [TSl or [TH) is give up forcibly. Figure [3] 
shows the process to combine a reasonably. It shows 
ai, a2 and 03, the combined a (acom.) and analytical 
a {(XAna.) of SAFl at z = 0, z = 1 and z = 2. It is found 
that there are evident calculation errors for ai, a2 and 
a^. The correlation coefhcients between ai, 2, 3, Com. and 
otAna. are given in table 2. However acom. is improved 
evidently, acom. match very well with a Ana., the cor- 
relation coefficients between acom. and aAna. are 0.99. 
Therefore the process to combine a is valid for SAFl 
(note this validation is not suitable for all cases of ex- 
trapolated field, but it can be valid on the whole, which 
will present next section). 

Figure S] shows the images of ai, 2. 3, Com. deduced 
from SAFl, AVI, BIE and Opt. fields at z = 1 and 
z = 2. There are evident deviations of a between ex- 
trapolated field and semi-analytical field. The correla- 
tion coefficients of ai,2,3,Com. between SAFl and each 
extrapolated fields are given in table 3. This low re- 
lation of as between extrapolated field and SAFl is 
caused mainly by the deviations between extrapolated 
field and SAFl, since from the above analysis it is 
found that the correlation coefficients between acom. 
and aAna. are already greater than 99% for SAFl. 



RSD = 



N - 



1 

— ^(a, -(a))2 



(17) 



Table 2 The correlation coefficients between 01,2, 3, Con 
and aAna. for SAFl. 





"Com. 


ai 


Q!2 


Q!3 


z=0 


0.99 


0.97 


0.93 


0.90 


z=l 


0.99 


0.92 


0.89 


0.91 


z=2 


0.99 


0.97 


0.90 


0.93 



As for the NLFF field, force-free factor a should be 
a constant along one special field line. Hence a physical 
quantity: relative standard deviation (RSD), which can 
estimate the deviation of a along field lines is studied. 
RSD is defined by formula (fT7|). where N is the num- 
ber of points calculated along field line. Figure (SKA) 
shows RSD of aAna. along each field lines for SAFl. 
Because the values of RSD approach zero {{\RSD\) is 
only 0.0056 and the points is very concentrated, which 
can give the conclusion that the field lines can be traced 
and the deviation of aAna. of SAFl is negligible. It 
means that SAFl satisfies the force- free equations very 
well. 

In Figure [SI the possibility function (PDF) of RSD 
of ai,2,3,Com. along some selected field lines are plotted 
for SAFl, AVI, BIE and Opt. fields. The reason for 
RSD of SAFl are plotted is that they can give an esti- 
mation on calculation errors. It can be found that PDF 
of RSD of a is very narrow and {\RSD\) is 0.09-0.16 for 
SAFl, this value can be taken as the calculation error 
of {\RSD\). For extrapolated NLFF fields, there are ev- 
ident deviations of a along field lines since PDF of RSD 
of a is not concentrated and {\RSD\) of a along some 
selected field lines is 0.96-1.19, 0.63-1.07 and 0.43-0.72 
for AVI, BIE and Opt. field, respectively. Addition- 
ally, the validation of combined is not evident for some 
cases, for example BIE extrapolated field (its RSD of 
Q^Com. is not the smallest one). In Figure [7l RSD vs 
field line length are plotted for SAFl, AVI, BIE and 
Opt. fields in order to find whether or not RSD are 
depend on field line length, where field line length was 
indicated by the number of points calculated along each 
field line. It can be found that RSD dose basically not 
depend on field line length. 

3.2 SAF2 

Figure [8] shows the field lines of SAF2, AVI, BIE and 
Opt. fields. Same as Figure [T] it can be found that 



Table 3 The correlation coefficients of Ofi, 2, 3, Com. between 
SAFl and the corresponding extrapolated fields. 









AVI 


BIE 


Opt. 


acom 


(z = 


1) 


0.36 


0.49 


0.67 


acom 


(z = 


2) 


0.42 


0.41 


0.59 


a2{z = 


= 1) 




0.43 


0.59 


0.62 


Q!2(z = 


= 2) 




0.46 


0.46 


0.54 


ai{z : 


= 1) 




0.33 


0.39 


0.55 


ai{z : 


= 2) 




0.38 


0.41 


0.56 


"3(2 = 


-1) 




0.65 


0.45 


0.53 


a3.{z : 


= 2) 




0.58 


0.52 


0.63 
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although there are some fine differences of field lines 
among these fields, the distributions of field lines of ex- 
trapolated fields basically match that of SAF2 at large 
scale. Like SAFl, we also calculate (F), F^ax and aj 
for each extrapolated field corresponding to SAF2. The 
results are given in table 4, here the pixel size is also 
assumed 1 arcsec. It can be found that for SAF2 angles 
between J and B are large than those for SAFl, even 
for semi-analytical field. The {F) and Fmax differences 
among extrapolated fields are negligible, however the 
aj of Opt. field is better than those of AVI and BIE 
fields. 

Same as Figure [2] for SAFl, Figure IH] shows the im- 
ages of ai, a2 and as, components of magnetic field 
and current for SAF2 aX z ~ 0. This also demonstrates 
that calculation errors are mainly located near where 
the magnetic field components reversed. Like Figure [3] 
for SAFl, Figure [To] shows ai, a2, as, acom. and cx.Ana. 
of SAF2 at z = 0, z = 1 and z = 2. Through calculation 
we get the correlation coefficients between ai^2, 3, Com. 
and a Ana., which arc shown in table 5. It shows that 
OiCom. evidently improve the computational precision 
of a for SAF2. 

Figure [TT] shows the images of a 1,2, 3. Com. deduced 
from SAF2, BIE, AVI and Opt. fields at z = 1 and 
z = 2. Comparing the results of SAFl, the consis- 
tency between as of extrapolated field and SAF2 is im- 
proved at some extent. The correlation coefficients of 
Q!i,2,3,Com. between SAFl and the corresponding ex- 
trapolated fields are given in table 6. This relative 
low relation of as between extrapolated field and SAF2 
is also caused by the deviations extrapolated field and 
semi-analytical field. 

Figure E^B) shows RSD of UAna. for SAF2, where 
{\RSD\) is 0.0023 and the points is also very concen- 
trated, which also indicate that deviation of aAna. of 
SAF2 is negligible. It means that SAF2 also satisfies 
force-free equations very well. 

Same as Figure El PDF of RSD of a along some 
selected field lines are plotted in Figure [T^] for SAF2, 
BIE, AVI and Opt. fields. From this figure, it can 
be found that calculation errors denoted by RSD, that 



Table 4 The results of the average {{F)) and maximum 
(Fmax) Lorentz-force and (Jj for SAF2 and the correspond- 
ing extrapolated fields. 





{F) 

{GyM)x 10-12 


^ max 

{GyM)x 10-" 




SAF2 


1.4 


13.6 


0.30 


AVI 


9.3 


25.7 


0.69 


BIE 


7.7 


32.4 


0.71 



Opt. 1.8 2.6 0.31 



deduced from SAF2, is smaller than those of SAFl. 
{\RSD\) is 0.13-0.11, 0.80-1.02, 0.67-1.34 and 0.33-0.55 
for SAF2, AVI, BIE and Opt. fields, respectively, which 
are better than those of SAFl. However the PDFs of 
RSD are also very wide for these extrapolated NLFF 
fields. RSD vs field line length for SAF2, AVI, BIE and 
Opt. fields arc plotted in Figure [T31 it can give the 
same results consistent with those of SAFl, that RSD 
do not depend on field line length either. 

4 Discussions and Conclusions 

In this paper, force-free factor a of NLFF field was 
mainly studied. The aim is to find, to what extent, a 
along a given field line can keep a constant, and to give 
an error estimation on a calculated from magnetic field. 

Through analysis, it is found that there are unavoid- 
able calculation errors for deducing a from magnetic 
field, the calculation errors are most likely to locate near 
where magnetic field components are reversed. ( | RSD \ ) 
of a along selected field lines is about 0.1 ^ 0.2 for semi- 
analytical fields, which can be considered as calculation 
errors of {\RSD\) caused by computation completely. 

It is found that there are obvious deviations on a of 
extrapolated fields from that of semi-analytical fields. 
The results of deviation of a along selected field lines 
arc as follows: For SAFl, {\RSD\) of a along selected 
field lines are about 0.96-1.19, 0.63-1.07 and 0.43-0.72, 
for AVI, BIE and Opt. fields, respectively. While for 
SAF2, they are about 0.13-0.11, 0.80-1.02, 0.67-1.34 
and 0.33-0.55 for AVI, BIE and Opt. fields, respec- 
tively. In both cases, it can be found that RSD of a 
along the selected field lines do not depend on field line 
length. 

Since RSD is a criterion for the deviation from a per- 
fect force-free state (the lower values indicate a more 
accurate force-free state) and RSD of Opt. field is 
less than those of AVI and BIE fields basically, so the 
performance of Opt. extrapolated fields is superior to 
other two extrapolated fields at some extent. In addi- 
tion, an interesting thing is that acom. is valid for all 
case of Opt. field, which is the same as semi-analytical 
fields. Previous studies mostly paid close attentions to 
the global or point to point properties for extrapolated 

Table 5 he correlation coefficients between 01,2,3, Com. 
and aAna. for SAF2. 





^Com. 


ai 


a2 


ag 


z=0 


0.96 


0.96 


0.92 


0.92 


z=l 


0.99 


0.94 


0.94 


0.94 



z=2 0.99 0.93 0.97 0.97 



fiedl, even if these properties are reasonable and accept- 
able, another property of force-free field (a should be 
an constant along field line) may not be satisfied well. 
For example, BIE method have given an constrain on 
the angles between B (calculated from Helmholtz equa- 
tion) and J (deduced from B) at each point, so the 
reasonability of global performance may be improved 
on the whole, but higher requirements (a should be 
an constant along field line) should be add for extrap- 
olation method. For the extrapolation errors of AVI 
method, there are two main effects, first is that it re- 
constructs the field by two field terms; second is that 
the singularity problems can not removal completely. In 
fact, the initial condition of potential field for Opt. field 
may give better results of RSD at some extent, however 
the superiority of Opt. field can also be studied from 
other aspects, such as the correlation coefficients of a 
between extrapolated field and semi-analytical field and 
other global properties. 
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Fig. 1 The magnetic field lines for each NLFF fields 
(SAFl, AVI, BIE and Opt. methods are labeled). 



Table 6 The correlation coefficients of 01,2,3, Com. between 
SAFl and the corresponding extrapolated fields. 





AVI 


BIE 


Opt. 


acomiz = 1) 


0.54 


0.77 


0.72 


aCom.(z = 2) 


0.59 


0.73 


0.68 


a2{z = 1) 


0.48 


0.71 


0.67 


a2{z = 2) 


0.54 


0.64 


0.54 


ai{z ~ 1) 


0.56 


0.65 


0.71 


ai(z = 2) 


0.51 


0.63 


0.66 


Q.z{z = 1) 


0.65 


0.58 


0.68 



az[z = 2) 0.58 0.55 0.65 









kit! 



Fig. 2 The images of a\, 02 and 03, the components of 
magnetic field and current [J^, Jy, J^) labeled in each frame 
for SAFl at 2 = 0. 
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Fig. 4 The images of acom. ,1,2,3 deduced from extrapo- 
lated fields at 2 = 1 and 2 = 2. (SAFl) 
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Fig. 5 RSD (Sa) of aAna. calculated from selected field 
lines for two semi-analytical fields (labeled). 
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Fig. 7 The scatter diagram of RSD {Sacam 123) ^^e 
length of field line calculated from some selected field lines 
for each NLFF field. (SAFf) 




Fig. 8 The magnetic field lines for each NLFF fields 
(SAF2, AVI, BIE and Opt. methods are labeled) 
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Fig. 11 The images of acom. ,1,2,3 deduced from extrapo- 
Fig. 10 The images of ai, ai and ai, oicom. and 01 Ana. at lated fields at 2: = 1 and z = 2. (SAF2) 
z = 0, z = 1 and z = 2 for SAF2. 
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Fig. 12 The PDF of RSD {Sacom. .1,2.3) calculated from 
some selected field lines for each NLFF field. (SAF2) 




Fig. 13 The scatter diagram of RSD ((5acom 1 3) 
length of field line calculated from some selected field lines 
for each NLFF field. (SAF2) 
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